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Abstract. We present results of several pointed X-ray broad band observations of the 'symbiotic X-ray binary' 4U 1954-1-31 
performed with the satellites BeppoSAX, EXOSAT, ROSAT, RXTE and Swift between October 1983 and April 2006. We also 
studied the RXTE ASM data over a period of more than 10 years, from January 1996 to October 2006. Light curves of all pointed 
observations show an erratic behaviour with sudden increases in the source emission on timescales variable from hundreds to 
thousands of seconds. There are no indications of changes in the source spectral hardness, with the possible exception of the 
RXTE pointed observation. Timing analysis does not reveal the presence of coherent pulsations or periodicities either in the 
pointed observations in the range from 2 ms to 2000 s or in the long-term RXTE ASM light curve on timescales from days 
to years. The 0.2-150 keV spectrum, obtained with BeppoSAX, is the widest for this source available to date in terms of 
spectral coverage and is well described by a model consisting of a lower-energy thermal component (hot diffuse gas) plus a 
higher-energy (Comptonization) emission, with the latter modified by a partially-covering cold absorber plus a warm (ionized) 
absorber. A blackbody modelization of our BeppoSAX low-energy data is ruled out. The presence of a complex absorber local 
to the source is also supported by the 0. 1-2 keV ROSAT spectrum. RXTE, EXOSAT and Swift X-ray spectroscopy is consistent 
with the above results, but indicates variations in the density and the ionization of the local absorber. A 6.5 keV iron emission 
line is possibly detected in the BeppoSAX and RXTE spectra. All this information suggests that the scenario that better describes 
4U 1 954-1-3 1 consists of a binary system in which a neutron star orbits in a highly inhomogeneus medium, accreting matter from 
a stellar wind coming from its optical companion, an M-type giant star. 
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accretion disks 



1. Introduction 

X-ray binaries are interacting systems composed of a com- 
pact object, neutron star (NS) or black hole (BH), which ac- 
cretes matter from a "normal" star orbiting at close distance. 
According to the mass of the donor star, these systems are clas- 
sified as High Mass X-Ray Binaries (HMXBs) or Low Mass 
X-ray Binaries (LMXBs). Among LMXBs, the mass transfer 
generally occurs through inner Lagrangian point matter over- 
flow from a Roche Lobe filling red dwarf (e.g.. White et al. 
1995; van Paradjs & McClintock 1995). 

However, there exists a very small subclass of LMXBs (3 
out of more than 150 objects: Liu et al. 2001; Masetti et al. 
2006a) that host an M-type giant mass-losing star. This red gi- 
ant is believed to feed the compact companion via its stellar 
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wind. By analogy with the symbiotic stars, which are systems 
composed of an M-type giant and a white dwarf (WD), these 
LMXBs are sometimes referred to as 'symbiotic X-ray bina- 
ries'. The protoype of this small subclass, the system GX 1h-4, 
is a well-studied case (e.g., Chakrabarty & Roche 1997); the 
second member, object 4U 1700H-24, is now receiving more 
attention (Masetti et al. 2002, 2006a; Galloway et al. 2002; 
Tiengo et al. 2005). 

The third member of this class, 4U 1954H-31, despite being 
relatively luminous in X-rays (it is an Uhuru source), has been 
neglected over the years. The object was first detected in the 
'70s in this spectral window by Uhuru (Forman et al. 1978) and 
subsequently hy Ariel V (Warwick et al. 198 1) in the framework 
of their all-sky surveys. Warwick et al. (1981) found the source 
at an average flux of Ri2x 10"^*^ erg cm^^ s^^ in the 2-10 keV 
band, with long-term variability of a factor of more than 20. 
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Table 1. Log of the X-ray pointings on 4U 1954+3 1 analyzed in this paper. 



Satellite 


Obs. start 
date 


Obs. start 
time (UT) 


Exposure 
(ks) 


On-source time (ks) 


EXOSAT 


22 Oct 1983 


21-50-00 


15.2 


14.7 


(GSPC) 






21:53:04 


14.8 








02 Jun 1985 


06:54:08 


16.3 


14.8 


(GSPC) 






06:57:20 


5.3 


5.2 


(ME) 






08:37:36 


10.2 


8.9 


(ME) 


ROSAT 


04 May 1993 


12:28:52 


98.1 


7.9 


(PSPC B) 


RXTE 


14 Dec 1997 


03:44:11 


15.0 


8.4 


(PCA) 






03:44:11 


15.0 


2.6 


(HEXTE) 


BeppoSAX 


04 May 1998 


22:05:14 


92.4 


20.0 


(LEGS) 






22:05:14 


92.4 


46.4 


(MEGS) 






22:05:14 


92.4 


20.0 


(HPGSPG) 






22:05:14 


92.4 


21.0 


(PDS) 


Swift 


19 Apr 2006 


21:02:10 


7.1 


2.5 


(XRT) 



*Results from this observation were partially presented by Gook et al. (1984) 



The first pointed observation was performed with EXOSAT 
in October 1983. The data (Cook et al. 1984) showed flaring 
behavior on 1-10 min timescales but no evidence of pulsa- 
tions or periodic variations; the 1-20 keV spectrum was mod- 
eled with either a powerlaw with photon index T ~ 1.5 or a 
bremsstrahlung model with kT ^ 36 keV. These character- 
istics led Cook et al. (1984) to conclude that 4U 1954+31 is 
probably an HMXB hosting a compact object accreting from 
a highly inhomogeneous wind from a companion star of early 
spectral type. 

Subsequently, 4U 1954+31 was observed in October 1987 
by Ginga in the 2-30 keV band (Tweedy et al. 1989). This ob- 
servation confirmed the timing characteristics found by Cook 
et al. (1984) and aflowed Tweedy et al. (1989) to model the 
spectrum of the source with a high-energy cutoff powerlaw 
plus a soft excess below 4 keV, heavily absorbed by interven- 
ing matter (possibly local to the source) partially covering the 
soft X-ray emission. These authors also found spectral vari- 
ability which they interpreted as due to changes in both the 
continuum shape and the absorption column. Tweedy et al. 
(1989) also (unsuccessfully) searched for the optical counter- 
part of this source, concluding that the M-type giant within the 
EXOSAT error circle is unlikely associated with 4U 1954+31 
and that the true counterpart should be an -emitting, highly 
reddened object (most likely a supergiant secondary hosted in 
an HMXB). 

Notwithstanding all of the above, which makes 4U 
1954+31 a potentially interesting target for X-ray studies, no 
refereed publications were devoted to this source, and no opti- 
cal counterpart has been associated with it, up to March 2006. 
After a hiatus of nearly 20 years, Masetti et al. (2006a), thanks 
to a pointed Chandra observation and a spectroscopic optical 
follow-up, showed that the M-type giant mentioned above is in- 
deed the optical counterpart of 4U 1954+31, and that this sys- 



tem is a symbiotic X-ray binary located at a distance d 
kpc. 

The work of Masetti et al. (2006a) shed light on this high- 
energy emitting source, and revived the interest in this object: 
a fast flare (lasting less than 2 hours) from 4U 1954+31 was 
detected with INTEGRAL in April 2006 and was communi- 
cated by Paltani et al. (2006); also, Corbet et al. (2006) found 
a 5-hour periodic modulation in the Swift BAT data acquired 
between December 2004 and September 2005. Mattana et al. 
(2006), using INTEGRAL observations and a subset of the data 
presented in this work, confirmed the NS spin period and the 
high absorption in the X-ray spectrum. 

In the present paper, we coUect aU of the unpubUshed 
archival X-ray observations on this system to fully characterize 
its high-energy behaviour and characteristics in both temporal 
(long- and short-term) and spectral scales. Particular attention 
is paid to the still unexplored parts of the X-ray spectrum of 4U 
1954+31, i.e. below 2 keV and above 30 keV, where the high 
sensitivity of BeppoSAX is of paramount importance to gauge 
the source emission. 

The aim of this paper is therefore to study the high energy 
properties of 4U 1 954+3 1 , and in particular the soft excess seen 
by Tweedy et al. (1989) below 4 keV. With this purpose, we 
have analyzed the data collected on 4U 1954+31 obtained over 
a time span of ~23 years through pointed observations per- 
formed with various satellites {EXOSAT, ROSAT, RXTE, Swift 
and BeppoSAX). We have also retrived the data collected with 
the RXTE AU-Sky Monitor (ASM), in order to have a better 
view of the source long-term X-ray aclivily. 

2. Observations and data reduction 

Here we describe the X-ray observations of 4U 1954+31 made 
with five different spacecraft. Table 1 reports the log of all the 
X-ray pointed observations presented in this paper. 
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2.1. BeppoSAX data 

This source was observed with the Narrow Field Instruments 
(NFIs) onboard BeppoSAX (Boella et al. 1997a) on 4 May 
1998. The NFIs included the Low-Energy Concentrator 
Spectrometer (LEGS, 0.1-10 keV; Parmar et al. 1997), two 
Medium-Energy Concentrator Spectrometers (MFCS, 1.5- 
10 keV; Boella et al. 1997b), a High Pressure Gas Scintillation 
Proportional Counter (HPGSPC, 4-120 keV; Manzo et al. 
1997), and the Phoswich Detection System (PDS, 15-300 keV; 
Frontera et al. 1997). LECS and MFCS were imaging instru- 
ments, whereas HPGSPC and PDS were non-imaging colli- 
mated detectors with fields of view of about 1?0 and 1?3, re- 
spectively. During all pointings the four NFIs worked nomi- 
nally and the source was detected by all of them. 

Good NFI data were selected from intervals outside the 
South Atlantic Geomagnetic Anomaly when the elevation an- 
gle above the earth limb was > 5°, when the instrument func- 
tioning was nominal and, for LECS events, during space- 
craft night time. The SAXDAS 2.0.0 data analysis package 
(Lammers 1997) was used for the extraction and the processing 
of LECS, MFCS and HPGSPC data. The PDS data reduction 
was instead performed using XAS version 2.1 (Chiappetti & 
Dal Fiume 1997). LECS and MECS data were reduced using 
an extraction radius of 4', centered on the source position; be- 
fore extraction, data from the two MECS units were merged. 

Background subtraction for the two imaging instruments 
was performed using standard library files, while the back- 
ground for the PDS data was evaluated from two fields ob- 
served during off-source pointing intervals. Due to the pres- 
ence of the strong X-ray source Cyg X-1 in one of the off- 
source pointings, only the other one was used for the back- 
ground subtraction. For the same reason, the background for 
the HPGSPC data was computed through an Earth-occultation 
technique (Manzo et al. 1997). 

Inspection of the MECS data shows the presence of a field 
X-ray source, located at ^21' from 4U 1954H-31, which we 
identify with the dwarf nova FY Cyg. By studying its X-ray 
spectrum (see Appendix) and extrapolating it to the ranges 
covered by the non-imaging instruments HPGSPC and PDS, 
we could determine the degree of contamination of the X- 
ray emission of this source on the hard X-ray spectrum of 4U 
1954H-3 1 . It is found that the X-ray emission of EY Cyg is less 
than 0.1% of the flux detected from 4U 1954H-31 above 7 keV 
(see Sect. 3.2), thus negligible. No contamination from EY Cyg 
is present in the MECS data because it appears fully separated 
from 4U 1954H-31 in the observation acquired by this imaging 
instrument. EY Cyg was not observed by the LECS as it was 
outside its field of view. 

We also checked for the presence of other contaminating 
X-ray sources in the PDS and HPGSPC fields of view. Only 
three RO SAT somces (Voges et al. 1999, 2000) were found at 
the edge of the PDS field of view, at distances of more than 
35' from 4U 1954H-3 1 . Assuming conservatively a hard spectral 
photon index F = 1 for all of them, we determined their total 
flux convolved with the PDS triangular response. We found that 
their contamination in the 15-150 keV band is less than 0.9% 



of the flux of 4U 1954H-31 in this band (see Sect. 3.2), and 
therefore negligible as well. 

In the same way, we evaluated the contamination of field 
X-ray sources on the PDS background measurement: 5 ROSAT 
sources were found in the off-source field used for background 
evaluation; it is found that these contribute less than 0.35% of 
the total measured background. Thus, again in this case, we can 
neglect this contribution. 

2.2. Other public X-ray data 

In order to present a thorough long-term analysis of the X- 
ray behaviour of 4U 1954H-31, we also retrieved from the 
HEASARC archiv43 all the unpublished X-ray observations 
of this source. These include (see Table 1) ROSAT, RXTE, 
EXOSAT and Swift data; we also re-examine the EXOSAT 
pointing already partially studied by Cook et al. (1984). 

Two pointed observations were performed on 4U 1954H-31 
with the EXOSAT ssLtellite (White & Peacock 1988). The first 
one (EXOSAT I) was acquired on 22 October 1983 while 
the second one (EXOSAT II) on 2 June 1985. Both observa- 
tions were made using the Medium Energy (ME; Turner et al. 
1981) proportional counter, in the 1-50 keV band, and with the 
Gas Scintillation Proportional Counter (GSPC; Peacock et al. 
1981), in the range 2-20 keV. However, due to their very low 
signal-to-noise ratio (S/N), GSPC data of EXOSAT II were not 
used in this analysis. 

The ROSAT satellite (Triimper 1982) was pointed at 
4U 1954H-31 on 4 May 1993 with the Position Sensitive 
Proportional Counter (PSPC; Pfeffermann & Briel 1986) unit 
B, sensitive in the 0.1-2.4 keV range. 

The RXTE sateUite (Bradt et al. 1993) observed the 
source on 14 December 1997; this satellite carries a 5-unit 
Proportional Counter Array (PCA; Jahoda et al. 1996), which 
is sensitive in the 2-60 keV energy range and allows a time res- 
olution of 1 /is, and a High Energy X-ray Timing Experiment 
(HEXTE; Rothschild et al. 1998) composed of two clusters of 
4 phoswich scintillation detectors working in the 15-250 keV 
band. 

This satellite also carries an ASV(1 (Levine et al. 1996) 
which regularly scans the X-ray sky in the 1.5-12 keV range 
with a daily sensitivity of 5 mCrab. Figure 1 reports the com- 
plete 1.5-12 keV ASM fight curve of 4U 1954H-31 starting on 
January 1996 up to October 2006, along with the times of the 
BeppoSAX, RXTE and Swift pointed observations. In order to 
clearly display the long-term trend in the X-ray emission from 
4U 1954H-31, each point in Fig. 1 is computed as the average 
of 5 subsequent measurements; given that the single original 
ASM points we retrived were acquired on a daily basis (i.e. 
they are one-day averaged measurements), the plot illustrated 
in Fig. 1 corresponds to a 5-day averaged X-ray light curve. 

Further X-ray data of 4U 1954H-31 were collected on 19 
April 2006 with the X-Ray Telescope (XRT, 0.2-10 keV; 

' available at: http : / /heasarc . gsf c . nasa . gov/ 
cgi-bin/W3Browse/w3browse . pi 

" ASM light curves are available at: 
http : //xte .mit . edu/ASM_lc . html 
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Fig. 1. 1.5-12 keV 5-day averaged RXTE ASM light curve of 4U 1954+31. One ASM count s^^ roughly corresponds to 13 
mCrab assuming a Crab-like spectrum. In the plot the times of the pointed RXTE, BeppoSAX and Swift are indicated by the 
vertical dashed lines (different dashes correspond to different spacecraft). A series of aperiodic increases in the X-ray activity 
from the source are noticed occurring on timescales variable from 200 to 400 days. 



Burrows et al. 2006) on board Swift (Gehrels et al. 2004). The 
data reduction was performed using the XRTDAS v2.4 stan- 
dard data pipeline package (xrtpipeline vO.10.3), in order 
to produce the final cleaned event files. 

As the XRT count rate of the source was high enough to 
cause pile-up in the data, we extracted the source events in an 
annulus, of 12" inner radius and 45" outer radius, centered on 
the source. The size of the inner region was determined fol- 
lowing the procedure described in Romano et al. (2006). The 
source background was measured within a circle with radius 
80" located far from the source. The ancillary response file was 
generated with the task xrtmkarf (vO.5.2) within FTOOLSH 
(Blackburn 1995), and accounts for both extraction region and 
PSF pile-up correction. We used the latest spectral redistribu- 
tion matrices in the Calibration Databas^ (CALDB 2.3) main- 
tained by HEASARC. 

We determined the source position using the 
xrtcentroid (vO.2.7) task. The correction for the mis- 
alignment between the telescope and the satellite optical axis 
was taken into account (see Moretti et al. 2006 for details). 
The Swift XRT position we obtained for 4U 1954+31 is the 
following (J2000): RA = 19^ 55'" 42!28; Dec = +32° 05' 45"5 
(with a 90% confidence level error of 3"6 on both coordinates). 

These coordinates are consistent with the (more accurate) 
ones measured in X-rays with Chandra and with those of the 
optical counterpart (Masetti et al. 2006a). We remark that the 



available at: 

http://heasarc.gsfc.nasa. gov/ f tools/ 
^ available at: |http : / /he as arc .gsfc.nasa. gov/ docs/hea 



uncertainty in the source coordinates as determined with the 
Swift XRT data may be underestimated because, as said above, 
the XRT observation is affected by pile-up and the PSF of the 
source image is thereby distorted. 

3. Results 

3. 1. Light curves and timing analysis 

The iJXr£ long-term (January 1996/October2006) 1.5-12 keV 
ASM light curve (see Fig. 1) of 4U 1954+31 shows an erratic 
behaviour of the source with sudden increases in its activity. 
The timing analysis of the 5-day averaged data has not revealed 
a coherent periodicity but only a variability on timescales be- 
tween 200 and 400 days. A similar indication is obtained using 
the dwell-by-dwell and the daily-averaged ASM data sets. 

We also studied the short-term X-ray light curves obtained 
with the pointed observations on 4U 1954+31 (see Fig. 2). The 
light curves have different bin times (from 50 s to 300 s) de- 
pending on their S/N. 

All light curves show an erratic flaring activity on 
timescales variable between 100 s and 1000 s, the longer ones 
corresponding to the more intense flares. In particular, EXOSAT 
I shows the most intense flaring activity among the observa- 
tions presented here, while during EXOSAT II the source was 
more active in the first part of the observation. 

As Cook et al. (1984) and Tweedy et al. (1989) pointed 
out some variability in the source hardness ratio as a function 
of the X-ray count rate, we investigated the presence of this 
^S£ibteadQttiglaM0itiltoi:i4tLI:A954iiai I . However, we did not find 
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Fig. 2. Background-subtracted X-ray light curves of 4U 1954+31 as seen during the pointed observations reported in the text. 
(Upper panels) 2-10 keV light curve of BeppoSAXMECS (left) and 0.1-2 keV light curve ofROSATVSVC B (right). The former 
is binned at 200 s, the latter at 300 s. (Central panels) 2-9 keV RXTE PCA (left) and 1-10 keV Swift XRT (right) light curves. 
Both Ught curves are binned at 50 s. (Lower panels) 0.8-8.9 keV light curves of EXOSATME observations I (left) and II (right). 
The former is biimed at 60 s, the latter at 100 s; for the sake of comparison of the activity of the source between the two EXOSAT 
observations, the latter panels were plotted using the same y-axis scaling. Given that, during the ROSAT observation, good source 
events were recorded only at the beginning and at the end of the pointing, we plot here these two time intervals. In all figures, 
times are in seconds since the beginning of the observation as reported in Table 1. 
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any change in the spectral hardness with X-ray intensity in any 
observation. Only RXTE data suggest a sUght hardening with 
increased source intensity. 

Timing analysis on the X-ray light curves from the various 
spacecraft was performed with the FTOOLS task XRONOS, 
version 5.21, after having converted the event arrival times to 
the solar system barycentric frame. The data were searched for 
coherent pulsations or periodicities but no evidence for them 
was found in the range 2 ms to 2000 s, although above 10 s any 
coherent behaviour could be significantly masked by the flaring 
activity of 4U 1954+31 (see Fig. 2). The power spectral den- 
sities obtained for each pointing are characterized by red noise 
and show no significant deviation from the l//-type distribu- 
tion (where / is the time frequency) typical of the 'shot-noise' 
behaviour. 

Likewise, a Fast Fourier Transform (FFT) analysis was 
made to search for fast periodicities (pulsations or QPOs) in 
the 1-2000 Hz range within the RXTE PCA data. During the 
observation, data in "Good Xenon" mode with a 2"^" s (i.e. 1 
yLts) time resolution and 256 energy bands were available. To 
search for fast pulsations and/or kiloHertz QPOs we made sev- 
eral FFTs of 16 s long data segments in the 2-10 keV and 10- 
60 keV energy intervals and with a Nyquist frequency of 4096 
Hz. We then calculated the statistics (Buccheri et al. 1983) 
on each ^^3000 s long subinterval of the RXTE PCA observa- 
tion (see Fig. 2, central left panel) and finally added together 
the results. No presence of QPO peaks or of coherent pulsa- 
tions was detected. 

3.2. Spectra 

In order to perform spectral analysis, the pulse-height spectra 
from the detectors of all spacecraft were rebinned to oversam- 
ple by a factor 3 the full width at half maximum of the energy 
resolution and to have a minimum of 20 counts per bin, such 
that the statistics could reliably be used. For all detectors, 
data were then selected, when a sufficient number of counts 
were obtained, in the energy ranges where the instrument re- 
sponses are well determined. In all cases we considered the 
average spectra of 4U 1954+31 given that no substantial vari- 
ations in the spectral shape during each pointed observation 
were suggested by the inspection of the color-intensity diagram 
of the source (See Sect. 3.1). 

Spectral analysis was performed with the package XSPEC 
(Arnaud 1 996) v 1 1 . 3 . 1 . In the broad band BeppoSAX fits, nor- 
maUzation factors were appUed to LEGS, HPGSPG and PDS 
spectra following the cross-calibration tests between these in- 
struments and the MEGS: these factors were fixed to the best 
estimated values (Fiore et al. 1999). Similarly, a constant was 
introduced between HEXTE and PGA in the RXTE spectra as 
well as between ME and GSPG in the EXOSAT spectra to take 
the different instrumental sensitivities into account. The spectra 
from the two HEXTE clusters were fitted independently and the 
normalization between them was left free, because of the pres- 
ence of a small systematic difference between their responses. 

To all spectral models tested in this paper (the best-fit mod- 
els for each satellite pointing are reported in Table 2), we ap- 



plied a photoelectric absorption column, modeled using the 
cross sections of Morrison & McGammon (1983; wabs in 
XSPEC notation), to describe the line-of-sight Galactic neutral 
hydrogen absorption towards 4U 1954+31. In the following, 
the acronym 'dof ' means 'degrees of freedom' . Table 3 reports 
values or 90% confidence level upper limits of the Equivalent 
Width (EW) of any Fe emission line at ~6.5 keV in the spec- 
trum of each pointed observation. 

We first considered the 0.2-150 keV BeppoSAX spectrum 
as it is the one with the widest spectral coverage. The initial at- 
tempt to fit this spectrum with a simple model such as a power 
law was unsatisfactory and indicated the presence of a soft ex- 
cess below ~2 keV and of a spectral turnover above 20 keV. 
The use of single-component models therefore proved to be in- 
adequate, giving unacceptable reduced values (>3). 

Mattana et al. (2006) described the X-ray spectrum of this 
source above 18 keV with a phenomenological model made of 
a cutoff power law. Our 15-150 keV PDS spectrum can be well 
fit with the same model and our results are compatible with 
theirs. When including the lower part of their broadband X-ray 
data, Mattana et al. (2006) used an absorbed cutoff power law 
plus blackbody model; however, they note that their spectral 
description is poor at low energies. We reached the same con- 
clusion when applying this same model to our BeppoSAX spec- 
trum, which means that a different and more physical modeUng 
should be used. 

We therefore assumed, as already applied to the X-ray 
spectra of this type of source (Masetti et al. 2002), a compos- 
ite model made of a Gomptonization component (COMPTT in 
XSPEC; Titarchuk 1994) and a thermal component to describe 
the emissions above and below 2 keV, respectively. We however 
found that, in order to describe the BeppoSAX spectrum at low 
energy, a complex absorption pattern was needed, with differ- 
ent absorption colunons acting on the two components. In par- 
ticular we found that, while on the thermal component a simple 
neutral hydrogen absorbing column was needed, a mixture of a 
partial covering absorber and of an ionized absorber (as mod- 
eled by Balucinska-Ghurch & McGammon 1992) was found 
acting on the Gomptonization component. We also noted that 
the thermal component could be described with an optically 
thin hot diffuse gas (mekal; Mewe et al. 1985) at relatively 
low temperature (fcTMKKAL ~ 0.07 keV); the thermal compo- 
nent could not be fitted with a blackbody or a disk-blackbody 
emission. 

The spectral fit is slightly improved (Ax^ = 4 for 2 dof 
less) by the addition of a narrow (with width a fixed at 0.1 
keV) Fe emission line at ~6.5 keV. In order to estimate the 
significance of this line, we used the F-test statistics for the case 
of an additional component (e.g., Bevington 1969): we find that 
the Une has a chance improvement probabiUty of 1 x 10~^. The 
BeppoSAX spectrum and the corresponding best fit model are 
reported in Fig. 3, upper left panel. 

To further study the presence of the soft excess found in the 
BeppoSAX spectmm, we analyzed the 0.1-2 keV spectral data 
of the ROSAT observation. Due to the fact that the ROSAT spec- 
trum has a smaller energy coverage and suffers from a lower 
S/N than the BeppoSAX one, we had to constrain the absorption 
and Gomptonization component parameters to the BeppoSAX 
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Fig. 3. X-ray spectra of 4U 1954+3 1 obtained from the pointed observations described in the text. The fit residuals using the best- 
fit model reported in Table 2 are also associated with each spectrum. (Upper panels) Be/5/7o5AX LECS+MECS+HPGSPC+PDS 
(left) and ROSAT PSPC B (right) spectra. (Central panels) RXTE PCA+HEXTE (left) and Swift XRT (right) spectra. (Lower 
panels) spectra of EXOSAT observations I (ME+GSPC data; left) and II (ME data; right). For observation II, the spectra of the 
high and low states (first 5 ks and remaining part of the observation, respectively; see Fig. 2, upper right panel, for reference) are 
shown. 



best-fit values in order to study the optically thin plasma emis- 
sion with a reasonable degree of precision. We find that the 
results (Fig. 3, upper right panel) are consistent with those of 
BeppoSAX. It should be noted that an acceptable description of 



the spectrum (with x^/dof = 17.7/18) could also be achieved 
without adding the partially covering absorption column; the 
parameters describing the low-energy thermal emission are in 
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Table 2. Best-fit parameters for the X-ray spectra of 4U 1954+31 from the pointed observations described in this paper. In all 
models in which a Fe emission line is needed in the fit, its width is fixed at the value cr = 0.1 keV. Frozen parameters are written 
between square brackets. Luminosities are not corrected for the intervening absorptions and are computed assuming a distance d 
= 1.7 kpc (Masetti et al. 2006a); in the cases in which the X-ray data could not completely cover the X-ray interval of interest 
for the luminosity determination, an extrapolation of the best-fit model is applied. Errors are at 90% confidence level for a single 
parameter of interest. Observations are reported in chronological order from left to right. 



model: WABS*[PCFABS*WNDABS(COMPTT + GAUSS) + MEKAL] 



Parameters'* 


EXOSATl 


EXOSAT 11 


ROSAT 


RXTE 


BeppoSAX 


Swift 




(1-15 keV) 


(1-15 keV) 


(0.1-2keV) 


(2.5-100 keV) 


(0.2-150 keV) 


(1.8-8.5 keV) 


X^/dof 


217.7/204 


86.0/85 


20.0/18 


112.6/105 


280.1/247 


53.4/44 


^WABS 


[1.5] 


[1.5] 


[1.5] 


[1.5] 


1 5+0.6 


[1.5] 


^PCFABS 






[25.4] 




25.4±1.6 


1 4+2.6 


CpCFABS (%) 






[95] 




95±2 


[95] 


jyWNDABS 






[3.0] 


4 1+0.2 






-BwNDABS (keV) 






[1.86] 


[1.86] 


1.86iri°3 




kTo (keV) 


0.56±0.04 


94+" "5 


[1.28] 


1 156+0"^^ 

l.iJO_(, 


1.28±0.06 




kT^- (keV) 


^•°O-0.1() 




[9.9] 




9.9±0.8 


sir 


T 


9.4l«i 


3.21-^ 


[2.8] 




9 Q+0.3 


4+5 


JS^COMPTT (XlO~^) 


501^3 


9 1+3.6 
^•^-2.1 


[8.2] 


5 0+°" 

-'•"—2.4 


8.2±0.8 


g+120 


feTMEKAL(xlO-^keV) 










6.81?:^ 




J^MEKAL 










200+5400 




Ep^ (keV) 






N/A 


o.46_o.ii 


6.46_o 11 




/Pe (xlO"*phcm-^ s"^) 






N/A 


1 Q+O.B 
^•^-0.4 


2.5±1.0 





X-ray luminosities'": 














0.1-2 keV 






4.1x10^^ 




2.0x10^^ 




2-10 keV 


1.8x10^^ 


3.4x10^'' 




7.9x10^" 


5.2x10^'* 


5.2x10^* 


10-100 keV 








1.2x10^5 


1.8x10^5 





"Hydrogen column densities are in units of lO'^'^ cm^'^. 
''In units of erg s~'. Values are not corrected for absorption. 
N/A: not applicable because of lack of spectral coverage. 



this latter case practically the same as those of the model with 
the partial absorption and reported in Table 2. 

The analysis of RXTE, Swift and EXOSAT spectra 
can instead allow one to monitor the behaviour of the 
Comptonization component and of the absorption column as- 
sociated with it as a function of time and of X-ray luminos- 
ity of the source. Using the model that best fits the BeppoSAX 
data we thus tried to describe the spectra of 4U 1954+31 se- 
cured with these X-ray satelUtes. However, as they do not 
cover the low-energy extreme of the X-ray band covered with 
BeppoSAX (i.e., below ~l-2 keV), they are insensitive to the 
presence of the mekal and to the Galactic litie-of-sight hydro- 
gen column. Thus, we modeled these spectra with an absorbed 
Comptonization model alone, and we fixed the Galactic Nh 
associated with the WABS component to the best- fit BeppoSAX 
value. 

The RXTE PCA+HEXTE spectrum, reported in Fig. 3, cen- 
tral left panel, apparently does not need the presence of a par- 
tial convering absorption, and only a warm absorber is re- 
quired. However, the fit is not formally acceptable (with x^/dof 
= 170.0/107) and shows large residuals around 6.5 keV. We 
therefore added a narrow emission line similar to that of the 
BeppoSAX best-lit model. This new description gives a Ax^ = 



57.4 for 2 dof less, meaning a chance improvement probabil- 
ity of 4xl0~^° according to the F-test statistics. The best-fit 
parameter values of the Comptonization component as deter- 
mined from the RXTE data are comparable to those measured 
from the BeppoSAX spectrum (see Table 2). 

EXOSAT pointings caught 4U 1954+31 at different emis- 
sion levels, with the source being on average six times brighter 
during observation I than in observation II. During both point- 
ings, the spectral description does not need any absorption be- 
sides the Galactic one. It is found that, while observation II 
data give Comptonization values broadly similar to those of 
the BeppoSAX observation (which was performed at a similar 
2-10 keV source flux), observation I shows that 4U 1954+31 
had a larger Compton cloud optical depth r and smaller elec- 
tron cloud and seed photon temperatures kT^- and kTo. We 
also find that the addition of an emission line at 6.5 keV is not 
needed in the best-fit model for both observations, with EW 
upper limits slightly tighter than, but broadly consistent with, 
those obtained by Gottwald et al. (1995) from the same data. 
The spectra of the EXOSAT pointings are reported in the lower 
panels of Fig. 3. 

Swift XRT data (Fig. 3, central right panel) span a smaller 
spectral range and have a lower S/N with respect to the 
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BeppoSAX, RXTE and fiXOSAT observations; nevertheless, the 
Comptonization model fits the Swift spectrum reasonably well, 
with parameter values (although affected by large errors) mid- 
way between those of the BeppoSAX and EXOSAT I best fits. 
The best-fit Swift spectrum does not need the presence of either 
an ionized absorber, or an iron emission line at 6.5 keV, and 
has a lower partially covering absorption coluimi density with 
respect to the BeppoSAX observation. 

We again explored the high-energy part of the BeppoSAX 
spectrum to look for the presence of a Cyclotron Resonant 
Feature (CRF) in absorption, which is considered one of the 
signatures of a highly magnetic NS hosted in an X-ray bi- 
nary system (e.g., Orlandini & Dal Fiume 2001 and references 
therein). We used a CYCLABS multiplicative model (Mihara et 
al. 1990; Makishima et al. 1990), but it did not statistically im- 
prove the best-fit model of Table 2. Among the other point- 
ings, thanks to its spectral coverage, only RXTE potentially al- 
lows one to perform a similar study. However, the lower S/N of 
HEXTE data did not allow us to explore the RXTE spectra. 

As a further test of the presence of a CRF, we also per- 
formed a Crab ratio of the BeppoSAX PDS spectrum. This 
technique was used to successfully pinpoint CRFs in numer- 
ous X-ray binary pulsars (see, e.g., Orlandini et al. 1998), and 
consists of obtaining the ratio between the count rate spec- 
trum of the source and that of the featureless power law spec- 
trum of the Crab. This ratio has the advantage of minimizing 
the effects due to detector response, non-perfect modeling of 
the spectral continuum and calibration uncertainties. Any real 
feature present in the source spectrum can then be enhanced 
by multiplying the ratio by the functional form of the Crab 
spectrum (a power law with photon index F = 2.1) and by 
dividing it by the functional form of the continuum adopted 
to fit the broad-band source spectrum, obtaining the so-called 
Normalized Crab Ratio (NCR; see Orlandini et al. 1998 for de- 
tails). 

As a result, the NCR performed on the PDS spectrum of 
4U 1954H-31 farther excludes the presence of any statistically 
significant absorption feature in the 15-150 keV range. This 
is at odds with the estimate of the magnetic field (B 10^^ 
Gauss) of the NS hosted in this system made by Mattana et al. 
(2006), as it would produce a clear CRF around 25 keV. The 
90% confidence level upper limit of the depth of any CRF at 
this energy is <0.25 assuming a feature width of 10 keV. 

4. Discussion 

The broadband, long-term X-ray study of 4U 1954H-31 pre- 
sented here, and performed through the analysis of data col- 
lected by five satellites, gave us the possibiUty to thoroughly 
study the long- and short-term high-energy behaviours of this 
bright but neglected source. In particular, with the widest ever 
broadband data coverage (0.2-150 keV) of this source, secured 
with BeppoSAX and shown here for the first time, we could 
characterize the X-ray spectrum of this object and check the 
best-fit model thus obtained against the observations acquired 
with RXTE, ROSAT, EXOSAT and Swift. Here we will discuss 
the results of the previous section and will compare them with 



Table 3. List of the 6.5 keV Fe emission line EWs, assuming 
a fixed line width of 0.1 keV, for the various X-ray observa- 
tions covering that spectral energy. Errors and upper limits are 
at a 90% confidence level. Pointings are listed in chronological 
order from top to bottom. 



Observation ^^(eV) 



EXOSAT \ 


<56 


EXOSAT \l 


<73 


RXTE 


4g+13 


BeppoSAX 


51 ±20 


Swift 


<1100 



the high energy behaviour of its twin, the symbiotic X-ray bi- 
nary 4U 1700-H24. 

4.1. Broadband X-ray properties of 4U 1954+31 

The temporal analysis of the fight curves of the RXTE ASM and 
of the pointed observations confirms that the source has pro- 
nounced long- and short-term variability, but also that it does 
not show any periodicity on either small (seconds to tens of 
minutes) or large (days to years) timescales. In particular, we 
did not find evidence for the 5 -hour periodic signal detected 
by Corbet et al. (2006) and Mattana et al. (2006) in any of the 
satelUte pointings. However, as can be seen from Fig. 2, this 
can be due to several factors: (a) the relatively short duration 
of the observations (in general lasting less than one single peri- 
odicity cycle); (b) the sparse light curve sampUng; and (c) the 
presence of large flares throughout each pointing. 

The presence of shot-noise stochastic variability (first 
found by Cook et al. 1984) is instead confirmed for the X-ray 
emission of 4U 1954H-31, while no short-period hardness ratio 
variations are detected during each pointed observation. These 
"colorless" short-term fluctuations in the X-ray flux may be 
produced by instabilities in the accretion process, or by inho- 
mogeneities in the velocity and/or density of an accreted stellar 
wind (or by both; e.g., Kaper et al. 1993). Thus, the observed 
random variability points to an explanation for this X-ray activ- 
ity as due to an inhomogeneous accretion flow onto a compact 
object (e.g. van der Klis 1995). 

The lack of evidence for X-ray eclipses in 4U 1954H-31 
suggests that either this system is viewed at a low (or interme- 
diate) inclination angle, or that the orbital period is quite long 
(hundreds of days). In fact, a hint of variations of the order of 
200^00 days is found in the secular ASM X-ray light curve. 

Concerning the spectral properties of 4U 1954H-31, one 
can see that the best-fit spectral model is typical of accret- 
ing systems hosting a NS (e.g., Paizis et al. 2006 and ref- 
erences therein). We cannot exclude the presence of a BH 
in this system; however, the temperatures associated with 
the Comptonization component are those generally seen in 
LMXBs with an accreting NS. The detection of a 5-hour mod- 
ulation (Corbet et al. 2006), interpreted as the spin of the ac- 
creting object, supports the accreting NS scenario. 
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The presence of a WD in accretion is basically ruled out in 
terms of both X-ray luminosity and spectral shape (de Martino 
et al. 2004; Suleimanov et al. 2005; Barlow et al. 2006). The 
X-ray emission from 4U 1954+31 is also much harder and 
stronger than that expected from the corona of a late-type giant 
star, this also being generally 4-5 orders of magnitude lower 
(Hiinsch et al. 1998) than the one detected in the pointed obser- 
vations presented here. 

As the source X-ray luminosity increases, the emission 
from the object becomes harder. This is consistent with a model 
in which enhancement of accreting matter onto a compact ob- 
ject is responsible for the observed spectral shape and variabil- 
ity. In addition, the hardening of the spectrum with increasing 
luminosity seems to indicate accretion of matter with low spe- 
cific angular momentum (i.e., from a wind), as suggested by 
e.g. Smith et al. (2002) and Wu et al. (2002) to explain a simi- 
lar spectral behaviour observed in other X-ray binaries. 

Comparison of X-ray spectral properties of the source 
among the different satelhte pointings shows a trend in the 
Comptonization parameters. In particular, an increase of r and 
a decrease of kTo and k.T,,- is found when the 2-10 keV X-ray 
luminosity increases above 10^^ erg s~^. At this luminosity, the 
source possibly undergoes a spectral state change. 

We found indications of complex and variable absorption 
around this X-ray source. A hint of this was already quahta- 
tively suggested by Tweedy et al. (1989), although their nar- 
rower spectral coverage did not allow them to explore this 
behaviour in detail. Indeed, the upper Umit to the interstellar 
Galactic absorption along the source Une of sight is estimated 
as N§'''^ = 0.84X 10^^ cm'^ (Dickey & Lockman 1990); from 
the BeppoSAX data we find that the Une-of-sight hydrogen col- 
umn is at least about a factor of two higher. The presence of a 
complex, variable and stratified neutral plus ionized absorbing 
medium may be explained by a NS accreting from a wind. 

Complex absorbers are sometimes seen at energies less than 
3 keV and can be described with a "two-zone" model (Haberl et 
al. 1989), where X-rays pass in series through a highly ionized 
absorber closer to the accretor and then a cold absorber farther 
away; 4U 1954+31 thus seems to belong to this source typol- 
ogy. This spectral behaviour is not unexpected from a source 
embedded in a stellar wind. In this sense, the detection of a rela- 
tively faint and narrow iron emission line suggests the presence 
of dense, ionized and radially infalling material in the close 
vicinity of the X-ray source. 

Likewise, the nondetection of an accretion disk is explained 
by the low angular momentum of the accreted stellar wind, 
which then falls nearly directly onto the NS surface without 
forming a large accretion disk structure. If the NS is magne- 
tized as suggested by the periodicity detected by Corbet et al. 
(2006), the accretion flow in the vicinity of the compact object 
is likely to be funneled onto the NS magnetic polar caps, thus 
leaving little room for the formation of a disk structure, albeit 
small. 

Concerning the magnetic field strength of the NS hosted 
in 4U 1954+31 we note that the nondetection of a CRF from 
this source is not unexpected, as other long spin period NSs 
do not show unambiguous CRF signatures (e.g., 4U 0114+65; 
Masetti et al. 2006b). The absence of a detectable CRF also 



casts some doubt on the possibility that this system hosts a 
highly-magnetized (B ~ 10^^ Gauss) NS as suggested by 
Mattana et al. (2006). 

A further indication that the accreting matter is flowing 
onto the NS polar caps via magnetic field confinement comes 
from the estimate of the size ro of the region emitting the 
Comptonization soft X-ray seed photons. Following the pre- 
scription by in 't Zand et al. (1999) for the computation of 
ro, and using the best-fit spectral parameters reported in Table 
2, we obtain that ro ^ 0.6 - 1.6 km, depending on the 
Comptonization state of the source. This estimate suggests that 
the area emitting soft seed photons on the NS covers only a 
fraction of its surface and it is comparable with the size of the 
base of an accretion column, which is wO.l times the NS radius 
(e.g., Hickox et al. 2004). 

From the above spectral information, we propose a model 
for this system similar to that developed by Hickox et al. (2004) 
for HMXBs hosting a NS accreting from the stellar wind of 
the companion. These authors indeed found that a soft excess 
around 0. 1 keV is a ubiquitous feature in the X-ray spectrum 
of HMXBs hosting a pulsating NS. Although the nature of the 
secondary is different in the case of 4U 1954+31, the proposed 
accretion mode (via dense stellar wind) is the same, therefore 
producing a similar behavior in the X-ray emission properties. 

Hickox et al. (2004) suggest that, for sources with luminos- 
ity Lx ^ 10^*^ erg s"^ (see right panel of their Fig. 3), die NS 
is embedded in a diffuse cloud. The inner parts of this cloud 
may be ionized, absorb the blackbody radiation photons from 
the NS surface and upscatter part of them via Comptonization 
mechanism. The (much cooler) external parts of this cloud also 
emit thin thermal plasma radiation, by collisional energization 
or reprocessing of a fraction of the hard X-ray emission com- 
ing from the Comptonization region inside the cloud by the 
optically thin external layers. The cloud is also responsible for 
the overall two-zone variable absorption detected in X-rays. 

In the Ught of all of these findings from our X-ray moni- 
toring campaign, we conclude that 4U 1954+31 is a symbiotic 
X-ray binary in which a mildly magnetized NS is embedded in, 
and accreting from, a complex and variable stellar wind coming 
from the secondary star, an M-type giant. 

4.2. Comparison with the system 4U 1 700+24 

We now draw a comparison between 4U 1954+3 1 and the other 
symbiotic X-ray binary 4U 1700+24. 

In terms of time scale variations of the X-ray emission, 
both systems show fast shot-noise variations, indicating that 
they are both powered by accretion onto a compact object, most 
likely a NS although up to now no evidence of a variability 
associated with a spin period was found in 4U 1700+24 (see 
Masetti et al. 2002 for details and an interpretation thereof in 
terms of system geometry). Galloway et al. (2002) suggested 
that the spin nondetection may be due to the fact that the NS 
could be slowly rotating. If true, this would be a further simi- 
larity between the two systems. 

Regarding the long-term variations, the behaviour seen in 
the ASM fight curve of 4U 1954+31 is reminiscent of that 
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shown by 4U 1700+24 and is interpreted as a modulation linked 
to the orbital motion of the system components (Masetti et al. 
2002; Galloway et al. 2002). For 4U 1954+31, however, we 
did not find variations which are as periodic and stable as those 
found in 4U 1700+24 (although see Mattana et al. 2006 for the 
possible presence of a 385-day period in 4U 1954+3 1). 

Spectral analysis of these two sources caught at different 
emission states shows that the Comptonization emission pa- 
rameters of 4U 1954+31, when observed at low luminosity 
(few 10^^ erg s~^ in the 2-10 keV band), are comparable to 
those shown by 4U 1700+24 during the phases of enhanced 
X-ray emission (when it reaches a similar X-ray luminos- 
ity, i.e., ^10^'' erg s"^ Masetti et al. 2002). So it seems that 
the Comptonization emission mechanism in these two sources 
forms a continuum depending on the luminosity state, thus sug- 
gesting a comparable overall accretion mechanism for these 
two X-ray binaries. 

The relatively high luminosity of 4U 1954+31, when com- 
pared to 4U 1700+24 (which is on average at luminosities 
which are 1-2 orders of magnitude fainter; Masetti et al. 2002), 
may indicate an intense stellar wind activity or, more likely, 
that the NS orbits much closer to the M-type giant than in the 
4U 1700+24 case. This would also explain the greater absorp- 
tion in the accreting flow, as in this way the NS hosted in 4U 
1954+3 1 would move into a denser medium. 

4U 1954+31, differently from 4U 1700+24, does not show 
the presence of blackbody emission from the accreting NS in 
its X-ray spectrum. This is reasonably explained by the fact 
that the former source shows the presence of a complex (and 
possibly variable) absorption around the X-ray emitter, which 
either absorbs the direct NS surface emission (at lower X-ray 
luminosities, when the accretion flow appears to be denser) or 
upscatters it via Comptonization processes (during the higher 
X-ray luminosity states, when the Compton cloud has a greater 
optical depth r). 

The two systems may however show a similar behaviour 
concerning the soft X-ray excess detected in both of them. 
Indeed, Tiengo et al. (2005) found in the XMM-Newton spec- 
trum of 4U 1700+24 an excess which they modeled with a wide 
Gaussian fine centered on 0.5 keV. A similar feature is present 
in the BeppoSAX spectrum of 4U 1954+31 presented here, and 
which we modeled as an optically thin plasma emission. It may 
thus be possible that the soft excess seen in 4U 1700+24 can be 
modeled in the same way. 

The absence of emission lines (or more generally of fea- 
tures connected with the accretion process) in the spectrum of 
tile optical counterpart of 4U 1954+31 (Masetti et al. 2006a), as 
in the case of 4U 1700+24, is naturally explained by the com- 
parison between the X-ray luminosity of 4U 1954+31 (f»10^^ 
erg s~^; see Table 2) and the total luminosity of its M-type gi- 
ant companion, -900 Lq (Lang 1992), i.e. 3.4x 10^^ erg s"\ 
for a M4 III star (Masetti et al. 2006a). As most of the lumi- 
nosity of the M-type giant is emitted in the optical and near- 
infrared bands, it can be understood that X-ray irradiation has 
httle effect on the secondary star. 

We speculate that the 4U 1954+31 system evolution is 
halfway between that of 4U 1700+24, where the compact ob- 
ject is accreting from a wind with lower density (probably be- 



cause it is orbiting farther away from the secondary), and GX 
1+4, where a high accretion rate onto the NS is apparent from 
X-ray and optical data (Chakrabarty & Roche 1997) implying 
that the latter system is tighter (with the secondary filhng its 
Roche lobe, thus producing an accretion disk around the com- 
pact object) tiian 4U 1700+24 and 4U 1954+31. 

5. Conclusions 

We have presented a broadband multi-spacecraft analysis of 
the symbiotic X-ray binary 4U 1954+31. This allowed us, 
for the first time, to explore in detail the X-ray characteris- 
tics of this 'forgotten' source. From the timing analysis we 
found that, while random X-ray variabiUty is present on long 
and short time scales, no periodic modulation is found. X- 
ray spectroscopy shows that the high-energy photons radiated 
from this source are produced by Comptonization plus thin 
thermal plasma emission. These two components, in particular 
the Comptonization, suffer from complex variable absorption, 
likely originating in the accretion fiow itself. An iron emission 
line is also possibly detected. 

From the interpretation of the above data and from the com- 
parison with similar sources, the scenario that we envisage is 
that of a mildly magnetized NS orbiting around an M-type gi- 
ant and accreting matter from the companion star stellar wind, 
which is inhomogeneous in terms of density and ionization de- 
gree, with the inner parts of the radial accretion flow possibly 
ionized by the NS high-energy radiation. We thus conclude that 
4U 1954+31 is a symbiotic LMXB which has several of the 
characteristics, in terms of accretion structure, usuaUy found in 
HMXBs. This likely led earlier studies to the conclusion (from 
the X-ray data analysis alone) that this source is an HMXB; 
only recently, thanks to combined optical and X-ray obser- 
vations, was it found that 4U 1954+31 is actually an LMXB 
(Masetti et al. 2006a). 

The soft excess observed in the X-ray spectrum of 4U 
1954+31 may be in the form of a group of discrete X-ray 
emission lines produced by light metals, rather than a con- 
tinuum emission. This is suggested by the fact that this ex- 
cess is modeled with an opticaUy thin thermal plasma model. 
Unfortunately, the relatively low spectral resolution of the data 
considered in this paper does not allow us to better explore this 
issue. Detailed modeling is needed to understand the nature of 
the observed complex spectral shape below 2 keV. To this aim, 
high-resolution spectroscopy with XMM-Newton or Chandra 
is, at present, the optimal choice. 
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Appendix A: The X-ray emission from thie fieid 
source EY Cyg 

As explained in Sect. 2.1, the Cataclysmic Variable EY Cyg 
is present in the MECS field as a faint X-ray source. In or- 
der to check for any possible contamination from this object 
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in the HPGSPC and PDS data, we extracted and analyzed the 
available information on EY Cyg as observed by the BeppoSAX 
MECS. Data were extracted from a region of 3' centered on the 
source; their reduction was similar to that appUed in Sect. 2. 1 
to the 4U 1954+3 1 MECS data, with the difference that, for the 
spectral analysis, we used the appropriate ancillary response 
function and background data for off-axis observations. 

The 2-10 keV light curve during the entire BeppoSAX ob- 
servation showed that EY Cyg had a steady behaviour, with a 
possible flare at half of the observation and lasting less than 
one hour (see Fig. A.l). The 5-10 keV / 2-5 keV hardness ra- 
tio of the source does not show significant variations during 
the BeppoSAX pointing: we therefore considered the 2-10 keV 
spectrum averaged over the entire observation. 
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Fig. A.l. MECS 2-10 keV countrate light curve of the source 
EYCyg. 



Table A.l. Comparison of the X-ray fluxes (expressed in erg 
cm-2 s-i) of EY Cyg and 4U 1954+31 in flie 7-30 keV and 
in the 15-150 keV bands, together with their percentage ratio, 
as measured during the BeppoSAX pointing described in Sect. 
2.1. 



Spectral 


EYCyg 


4U 1954+31 


Ratio 


range (keV) 






(%) 


7-30 


5.5xlO~" 


4.1x10-1° 


0.13 


15-150 


1.5x10 


3.8x10 


0.039 



Following the approach in Baskill et al. (2005), we fit the 
MECS spectral data using an absorbed optically thin thermal 
plasma (MEKAL model in XSPEC; Mewe et al. 1985). As no 
spectral information was available below 2 keV in the MECS 
data, we fixed the hydrogen column density to the value found 
by Baskill et al. (2005), i.e. 2.8x10^1 cm-^. The averaged 
spectrum was thus satisfactorily fitted (x^/dof = 8.3/14) with 
a kT = 8^4^ (see Fig. A.2), in broad agreement with the re- 
sults of Raskin et al. (2005). The best-fit model implied 7-30 
keV and 15-150 keV fluxes for the EY Cyg emission in the 
HPGSPC and PDS ranges as shown in Table A.l. 

As one can see from this Table, the comparison between the 
fluxes of 4U 1954+31 and EY Cyg indicates that any contami- 
nation from this field source on the spectrum of 4U 1954+3 1 is 
negligible. 



EY Cyg 

MECS data 




2 5 
Energy (keV) 



Fig. A.2. MECS 2-10 keV spectrum of the source EY Cyg, 
averaged over the BeppoSAX observation. Data (in the upper 
panel) are well fitted using a single-temperature opticaUy thin 
thermal plasma model with kT ~ 8 keV. In the lower panel, 
the best-fit residuals are shown. 



